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Abstract 

Recently, the method of one-dimensional maps was introduced as a means of generating excep- 
tional discretisations of the ^-theories, i.e., discrete </> 4 -models which support kinks centred at a 
continuous range of positions relative to the lattice. In this paper, we employ this method to obtain 
exceptional discretisations of the sine-Gordon equation (i.e. exceptional Frenkel-Kontorova chains). 
We also use one-dimensional maps to construct a discrete sine-Gordon equation supporting kinks 
moving with arbitrary velocities without emitting radiation. 
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I. INTRODUCTION 



Discrete analogs of nonlinear evolution equations have been the subject of intense inves- 
tigation over the last 15 years. A great deal of insight has been gained into the properties 
of the discretised nonlinear Schrddinger, Landau-Lifschitz, Korteweg-de Vries, 4 - and other 
equations which were originally introduced in the context of continuous nonlinear media. 
As for the discrete sine-Gordon equation that we study in this paper, it preceded the ap- 
pearance of its continuum counterpart in the physics literature. The equation dates back to 
1938 when it was proposed by Yakov Frenkel and Tatyana Kontorova to model stationary 
and moving crystal dislocations [l]. 

The original Frenkel-Kontorova model consisted of a chain of harmonically coupled atoms 
in a spatially periodic potential: 



Here 9 n is the position of the n-th atom in the chain, 1/h 2 is a coupling constant, and 
the overdots indicate differentiation with respect to time: 9 = d 2 9/dt 2 . An alternative 
interpretation of Eq. (CQ) is that of a chain of torsionally coupled pendula, with 9 n being the 
angle the n-th pendulum makes with the vertical. Finally, Eq.([T]) can be seen simply as a 
discretisation of the sine-Gordon equation 



which was conceived for the numerical simulation of this partial differential equation. 

Since its original inception as the Frenkel-Kontorova model, the discrete sine-Gordon 
equation (CQ) has reappeared in a great number of physical contexts, including domain walls 
in ferro- and antiferromagnetic crystals, charge-density waves in solids, crowdions in metals, 
vortices in arrays of Josephson junctions, incommensurate structures in metals and insu- 
lators, and nonlinear excitations in hydrogen-bonded molecules. (See for review and 
references.) The equation has also been generalised in a variety of ways. In the present 
paper we study, systematically, two classes of such generalisations. In the models of the 
first class the main part of the intersite coupling is still harmonic, as in the original Frenkel- 
Kontorova model, but in addition there is an anharmonic part of interaction arising from the 
modified periodic potential. These types of models are of interest primarily in the station- 
ary case where they define nontrivial systems of statistical mechanics (systems with convex 
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(1) 



9 = 9 XX — sin#, 



(2) 
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interactions) [3j, |4( . The stationary discretisations in this class have the form 

T^(#n+i — 2# n + Q n -\) = f(O n -i, 9 n , n+ i), (3) 

where the function / (not necessarily a periodic function) reduces to sin 9 in the continuum 
limit: 

f(9n-i,@n,9n+i) — > sin # n as 9 n -l, 9 n+1 — > n . 
The other class consists of all-periodic discretisations: 

^ sin(0 n+ i - 29 n + n _!) = /(^i, 9 n , 9 n+1 ). (4) 

Here / is a periodic function of each of its three arguments, such that f(9 n , 9 n , 9 n ) = sin9 n . 
Equations of this type govern arrays of electric dipoles or magnetic spins in which the inter- 
actions between the neighbouring elements are characterised by the trigonometric functions 
of the corresponding angles. These models are commonly referred to as the "sine lattices" 
j^]. In the stationary case, examples of the sine lattices include the usual and the chiral one- 



dimensional XY model in the magnetic field 



BE 



6|. In the time-dependent setting, the sine 



attices were used to model the rotational dynamics of methyl groups in 4-methyl-pyridine 
3], CH 2 units in crystalline polyethylene [8] and bases in a DNA macromolecule Q, 

to study conformational defects in polymer crystals [lO] and nonlinear waves in chains of 

electric dipoles 



In addition to (J3j) and (Hj), we also consider some other discretisations which have the 
special property of exhibiting exact solutions. 

In most physical applications of the sine-Gordon theory, both continuum and discrete, 
the central role is played by its solitary-wave solution, called a kink. The kink represents 
a dislocation in the crystal, a 27r-twist wave in the chain of pendula, and a quantum of 
magnetic flux in a long Josephson junction. In the continuum model (J2J), the kink solution 
is available explicitly; in particular, the stationary kink has the expression 

9{x) = 4arctan [exp (x — x^)}. (5) 

The stationary kink (jHJ) depends on a single parameter, the position of its center x^°\ which 
can be varied continuously: — oo < x^ < oo. In generic discretisations of the sine-Gordon 



theory, however, the kink can only 



two neighbouring sites 
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)e centred on a lattice site or strictly midway between 



151 ] . Mathematically, this is a consequence of the breaking 



of the translation symmetry of the continuum model. The physical interpretation is that 
the discrete kink can only remain stationary when placed at a minimum or a maximum of 
the so-called Peierls-Nabarro barrier, a periodic potential induced by the discretisation of 



equation 



13 



Speight and Ward 17j were the first to realise that the breaking of the translation sym- 



metry does not necessarily preclude the existence of a one-parameter family of stationary 
discrete kinks with an arbitrary centring relative to the lattice. In other words, despite 
not being translation invariant, the lattice equation may support a kink solution which de- 
pends on a continuous translation parameter. This "spontaneous symmetry restoration" is 
a nongeneric phenomenon which may only occur in isolated, or exceptional, discretisations 
of the sine-Gordon model. Physically, it implies that the discretisation does not induce the 
Peierls-Nabarro barrier, or that the barrier is transparent to kinks. Flach, Zolotaryuk and 
Kladko have discovered a similar phenomenon in a class of discrete Klein-Gordon systems 



with nonlinearities of a special form 18 



The classification of the exceptional discrete sine-Gordon equations, i.e., equations sup- 
porting families of stationary kinks with a continuously variable position relative to the 
lattice, is of fundamental interest. Firstly, the discrete kinks tend to be more mobile in 
exceptional discrete models. There are some isolated velocities at which kinks in the ex- 
ceptional models may slide, i.e. travel without losing energy to radiation [19(. For other 



velocities, the moving 



in generic systems 



dnks do radiate but the amplitude of radiation is much smaller than 



17 



19 



20 



21] . In addition, the collisions of kinks were reported to 



be more elastic in the exceptional models 22j. We also show in this paper that some ex- 
ceptional systems admit time-dependent versions which support sliding kinks with arbitrary 



velocities. Furthermore, there are indications 



23j that all exceptional discretisations possess 



a conservation law: they conserve either energy or momentum. Therefore, exceptional dis- 
crete models appear to be "better" approximations of the partial-differential equation (J2]) - 
at least as far as the kink solutions are concerned — as they preserve important properties 
of the continuum model. 

The objective of the present paper is to identify exceptional Frenkel-Kontorova models 
within the families dSj) and (jl]). We are not going to attempt a complete classification here; 
instead, we focus on identifying simple particular cases which may be of practical use in 
future. We also construct two discrete models with exact (stationary and moving) kink 
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solutions. 

An outline of the rest of the paper is as follows. In the next section (section [Til), we 
present the method of one-dimensional maps as applied to discrete sine-Gordon equations. 
In section IHIl the method is used to identify simple exceptional discretisations of the form 
([3]) involving ratios of trigonometric and linear functions. The symmetric maps found in this 
section have one further use; in section ITVl we utilise them to construct purely trigonometric 
discretisations [of the form (j3J)]. In the subsequent sections we present discrete sine-Gordon 
equations with exact stationary (section [V} and moving (section IVip kink solutions. Finally, 
several concluding remarks are made in section IVHI which summarises the results of this 
study. 



II. THE METHOD OF ONE-DIMENSIONAL MAPS 



To derive an exceptional discrete sine-Gordon model, Speight and Ward used the Bogo- 
molny energy- minimality argument 13]. The energy minimality requirement has naturally 
led them to consider a one-dimensional map rather than the original, second-order, dif- 
ference equation. In the follow-up work ]24], Speight utilised the energy-minimising map 
(the Bogomolny map) to prove the existence of a one-parameter family of kinks for their 
discretisation of the 4 -theory, i.e. the exceptionality of their 4 -model. 

A further insight was due to Kevrekidis |25l| . Inspired by Herbst and Ablowitz' results on 
the discrete nonlinear Schrodinger equation 26) , Kevrekidis reformulated the exceptionality 



of a stationary discrete Klein-Gordon equation as the existence of a two-point invariant. He 
also provided two phenomenological recipes of construction of stationary discretisations with 
such invariants. Thus, the existence of a two-point invariant replaced the energy minimality 
requirement as the crucial property of exceptional discretisations. 

The universality of one- dimensional maps as generators of translationally-invariant fami- 



lies of solutions has been fully realised in Ref. 



271 ] . Instead of trying to identify discretisations 



exhibiting a two-point invariant, the authors of [27j proposed to generate exceptional dis- 
cretisations departing from a postulated map. (That is, the two-point invariant has now 
become a starting point rather than the final objective of the analysis.) In this way, the 
classification of exceptional discretisations has been reduced to the classification of one- 
dimensional maps. This will be our approach in this paper as well. 
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We start by considering a discrete sine-Gordon equation of the form and assume that 
the corresponding stationary equation, 

J^(9 n +1 — 2#n + On-l) = f(0n-l, 6n, @n+l), (6) 

has a solution of the form 9 n = g(nh), where the continuous function g(x) is defined for 
— oo < x < oo and is monotonically growing, with g(—oo) = and g(oo) = 2tt. Since n does 
not appear in equation ([6]) explicitly, from the existence of the above solution it follows that 
equation (El) also has a whole family of solutions 9 n = g(nh — x^), with any real x^, and 
therefore, that the model ([6]) is exceptional. For each x^°\ the solution 9 n = g(nh — x^) 
represents a discrete kink; if we interpret values x n = nh as positions of the lattice sites on 
the x-axis, the kink 9 n appears centred on the point x = x^ . It is important to emphasise 
that we do not need to know an explicit form of g; all we need to know is that a function 
with these properties exists (for example, as an implicit function). 

As g(x) is a monotonically growing function, we can invert it to obtain 9 n+ \ = g{g~ 1 {9 n ) + 
h) = F(9 n ). Since g(x) is defined for all real x, the function F(9 n ) is defined for any 9 n . 
Thus the fact that the discretisation is exceptional implies that the kink solution 9 n satisfies 
a one- dimensional map [27|]. The opposite is also true. Namely, assume Eq.© results from 
the iteration of a one-dimensional map 9 n+ \ = F(9 n ) (in a similar way as a second-order 
differential equation can be derived by differentiating a first-order one). In addition, let 
the function F be such that F{6) > 9 for any 9 between and 2n, whereas F(0) = and 
F(2n) = 2n. A simple cobwebbing argument shows then that for any 9q within the range 
< 9 < 2ir, the map generates a discrete kink solution 9 , Oi, .... Therefore we have 

a one-parameter family of kinks and so Eq. © represents an exceptional discretisation. 

This observation implies that we can find exceptional discretisations of the sine-Gordon 
equation by considering a one- dimensional map 

9 n +i — 9 n = hH(9 n+ i, 9 n ), (7) 

where the function H satisfies several requirements. First of all, it should satisfy the condi- 
tion 

Q 

H(6, 9) = 2 sin - (8) 
2 

which ensures that the map ((7j) reduces to equation 

Q 

0a = 2 sin- (9) 
6 



in the continuum limit (where 6 n+ i — 9 n — > hO x ). Equation (P) is the Bogomolny equation for 
the stationary continuum sine-Gordon theory: the sine-Gordon equation 9 XX = sin9 follows 
from Eq.([9]) by differentiation, while its kink solution (jSJ) is simultaneously a solution of 
Eq.Q. Therefore, the condition (jSJ) selects maps which generate discretisations of the 
sine-Gordon rather than some other equation. Our second requirement is that H(6 n , 9 n +i) 
should be bounded and positive for all pairs of 9 n and 6 n+ \ with sufficiently small value of 
the difference \6 n+ i — 9 n \ (where < 9 n , 9 n+ i < 2ir). Using this property of H, assuming 
that h is sufficiently small, and invoking the implicit function theorem, we can show that 
Eq.Q defines, for any < 9 n < 2tc, a function 9 n +i = F(9 n ), with 9 n+ \ > 9 n . Thus Eq.((7j) 
will give rise to an exceptional discretisation of the sine-Gordon equation. 

This discretisation results from squaring both sides of (jZj) and subtracting the square of 
its back-iterated copy, 

9 n — 9 n -i — hH(9 n , 9 n _i). 



This yields [271 ] 



9 n +i — 29 n + 9 n -\ _ H 2 (9 n+ i, 9 n ) — H 2 (9 n , 9 n -\) 

If H is symmetric [i.e. invariant under the permutation of its arguments: H(x, y) = H(y, x)], 
the numerator in ( JTU1) vanishes whenever the denominator equals zero and hence the right- 
hand side of fllOp is nonsingular. Thus the classification of exceptional discretisations reduces 
to the classification of all symmetric functions H(x, y) with the above properties. The next 
section summarises results of this analysis. 

III. RATIONAL-TRIGONOMETRIC DISCRETISATIONS 

A. H 2 (x,y) = T(x) + T(y) 

The simplest posibility is to let H 2 (x,y) = J-{x) +J 7 (y). From H(x,x) = 2sin(x/2) it 
follows that T{x) = 2sin 2 (ic/2) and so 

H 2 (x,y) =2sin 2 ^ + 2sin 2 |. (11) 



This function gives rise to one of Kevrekidis' discretizations 251 ]: 



9 n+l - 29 n + fl w _! _ cosg n+ i - cosg n _i 

h 2 9n+l — On~l 
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or, equivalently, 

9 n +i — 26> n + On-i n sin [(9 n+ i — n _i)/2] . f 9 n+ i + 9 n -i , , 

S» = 2 tu.-<l-. sm l 2 (12b) 

B. H 2 (x,y) = F(x + y) 
Letting H 2 (x,y) = T{x + y) and substituting into the continuum limit condition, yields 

H 2 (x,y) = 4sin 2 ^±l. (13) 

This symmetric function generates the discretisation of the form 

9 n +i — 20 n + 0n-i ,sin [(0 n+ i — n _i)/4] . f 9 n+ i + 20 n + n _i\ , 
W = 4 fl^-fl^ Sm I 1 ) ■ (M) 

C. F 2 (x,y) =F{x)F(y) 

Another simple possibility is to assume that H 2 (x,y) = J r {x)J z {y). From the continuum 
limit we obtain T{x) = 2sin(x/2) and so 

H 2 (x, y) = 4 sin ^ sin — . (15) 

By substituting into equation ffTUl) we find the following discretisation: 

9 n +i — %9 n + 9 n -i sin[(0 n+ i — n _i)/4] . / n \ ( 9 n+ i + n -i\ . . 

sm — cos . (16) 



h 2 9 n+x -9 n - x \2J V 4 

D. J ff 2 (x,y) = [^(x)+^(y)] 2 

Considering the symmetric function of the form H 2 (x,y) = \!F(x) + jF(y)] 2 , we obtain 
from the continuum limit: 

H 2 (x,y) = (sin | + sin |) . (17) 
This gives rise to the following discretisation: 

9 n +i — %9 n + 9 n _ x ^ sin [(9 n+1 — g n _x)/4] ^ n+1 + n _ a 



/> 2 " 9 n+1 -9 n „ 1 V 4 



/ 9 n +\ . 9 n . 9 n _i . 
x (sin — + 2siiiy + sm— ]. (18) 
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A simple symmetric generalisation involving two functions of a single argument, say T{x) 
and Q(x), is H 2 (x,y) = + Q{x)Q{y). Setting x = y yields 

F\x)+g\x) = 4sin 2 |. 

One possibility here is to assume that the functions T and Q have the form 

T[x) = 2r](x) sin — , Q{x) = 2£(x) sin — , 

where rj and £ satisfy rf(x) + £ 2 (x) = 1. The simplest trigonometric choice for rj and £ is 

rj(x) = sm(ax), £(x) = cos(ax), 

where a is a parameter. Taking, for instance, a — |, gives us the following expression for 
H 2 (x,y): 

H 2 (6 n+1 , 9 n ) = 4 sin ^±1 sin | cos (19) 

The function ( fl9l) is obviously positive for \6 n+ \ — 9 n \ < it and therefore the resulting 
discretization, 

n +i — 26> n + dn-i .sin \{0 n+ i — 9 n -i)/2] . / 9 n \ ( 9 n+ i — 6 n + 0„_i\ 

~ 4 ~7i n sm I — 1 cos , (20) 



h 2 9 n+ i — 6 n _i \ 2 J \ 2 

is exceptional for sufficiently small h. 

Another simple symmetric combination of two functions of a single argument, is 
H 2 (x,y) = J 7 (x)Q(y) + !F(y)Q(x); however, this H 2 gives rise to the discretisation that 
we have already identified, Eq. (fT2l . 

F. More complex symmetric functions 

It is not difficult to construct more examples of symmetric functions H(x,y), with in- 
creasing complexity. One possibility is to take 

N 



H 2 (x,y) = ^2T n {x)F n {y), 



n=l 



where T n {x) (n = 1,2, ...,N) are appropriate trigonometric functions. Another symmetric 
combination is 

N N 

-E\n+1\ ) > 

n=l n=l 

where \n\ = n mod 2, and x\ = x, x^ = y. 



IV. PURELY TRIGONOMETRIC DISCRETISATIONS 



Our original one- dimensional map (J7|) can be modified to produce new, periodic, dis- 
cretizations. Instead of (J7J), we consider the map 

lsm 6n+1 ~ 9n =hH(9 n+1 ,9 n ), (21) 

where H(9 n+ i,9 n ) is a trigonometric function of its arguments and I is a positive integer. 
Subtracting from the square of (ED) the square of its back-iterated copy yields the discrete 
model 

P_ . 6 n +i — 20 n + Qn-1 _ H 2 (9 n+1 , 9 n ) — H 2 (9 n , fl n _i) 

-on - _ ______ . l ) 

As in Eq.(J7J), we assume that H(9 n+ i, 9 n ) is positive, symmetric and has the continuum limit 
H(9,9) = 2sin(0/2). For sufficiently small h and \9 n+ i — 9 n \, Eq. (j2ip defines an implicit 
function 9 n+ \ = F(9 n ), with 9 n+ \ > 9 n . Consequently, the discretisation (l22il is exceptional. 

The discretisations (1221) are different from those in (fTOl) in that every term in ( 1221) is 
periodic in each of its three arguments, 9 n _i, 9 n and The models of the form (122]) 

find their applications in the description of coupled chains of elements where each element is 
characterised by a periodic variable (an angle) and the coupling of elements does not violate 



this periodicity. One example is given by the Speight- Ward discretisation [13] 



4 . / 9 n+ \ — 29 n + 9 n -i \ . / 9 n+ i + 29 n + Q n -\ \ 
h? Sm { 1 ) = Sm { 1 ) ' (23) 

which has recently been shown to describe chains of electric dipoles constrained to rotate 



in the plane containing the chain Another example is the one-dimensional chiral XY 
model 

sin(x„+i - Xn -7) - sin(x„ - Xn-i ~ 7) = K sm(px n ) (p=l,2,...), (24a) 
or, equivalently, 

P . fOn+l - + &n-l \ ( ®n+l ~ #n-l \ . a fnA ,s 

sm cos 7 = sinfc)„, (24b) 



h 2 

where 9 n = pXn, i — 2p and h 2 = pK. The chiral XY model (j2*4*|) describes arrays of spins 
with the nearest-neighbour interactions in an external magnetic field {3, (|. It is used to 
model helimagnetic materials, discotic and ferroelectric smectic liquid crystals, crystalline 
polymers, thin magnetic films and Josephson junction arrays. 
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The classification of discretisations of the form (1221) reduces to the classification of all 
possible symmetric functions H(x, y) — the task completed in section ITTT1 above. Each func- 
tion H 2 (x,y) identified in section [TTT1 gives rise to a number of purely periodic discretisation 
of the form ( 1221) . with various £; that is, each rational-trigonometric exceptional model ( [TO]) 
has a set of purely trigonometric counterparts ( 1221) . We will restrict ourselves to the simplest 
representative(s) of these sets by choosing appropriate value(s) of £. The resulting models 
can be summarised as follows. 

Picking the symmetric function (TTTj) of the section [III Al and letting I = 2, gives rise to a 
very simple exceptional discretisation 

2 . [0 n +X — 2# n + 0n-l\ • {@n+l+0n-l\ /oc n, 

sin = sin . (25) 



h 2 \ 2 J \ 2 

If, instead, we took I = 4, we would obtain a slightly more complicated model: 

4 f 9 n+ i — 29 n + 9 n ^i\ . (9 n+1 + 9 n _i\ (9 n+ i — 9 n _ 1 \ 
- sin { i j = sm { j cos ^ i j . (26) 

Finally, if we "extend" the symmetric function ( [111 by adding a term that vanishes in the 
continuum limit, 

H 2 (x, y) = 2 sin 2 - + 2 sin 2 - + 2 sin 2 



2 ' "~ 2 ' 2 ' 

then, keeping i = 4, we will arrive at a (still reasonably simple) exceptional model 

- sm ^ i j = sm ^ j cos ^ i j cos -. (27) 

Next, choosing the symmetric function (fl3l) of the section Ull Bl and letting I = 4 yields 
Speight and Ward's model, Eq.([23l). On the other hand, taking the symmetric function ( fT5l) 
of section IIII CI and letting i = 4, gives the discretisation 

- sm ^ - J = 2 Sm U J C ° S I 4 J ■ m 

Eq. (f2"Hl) reduces to Eq. ([2Hl) with the lattice spacing constant h = h{l + /i 2 /4)~ 1//2 , if we use 
an identity 

2sm j cos ^ i j = sm ^ i j - sm ^ 

'(29) 

Picking the symmetric function of the form (1171) from section IIII D\ and letting £ = 4, 

gives an exceptional discretisation 

. , 9 n+ i — 29 n + # n _! \ 1 ( 9 n+ i + 9 n _i \ ( . 9 n+ i . 9 n . # n _i 
sm = - cos sm — h 2 sm — + sm 



h 2 \ 4 ) 2 \ 4 J V 2 2 2 

' (30) 
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Finally, choosing the symmetric function in the form (1191) from section IIIIEI and letting 

i = 2 produces an exceptional model 

2 . / n+ i — 29 n + 9 n -i\ . ( 9 n \ ( 6 n+ i — 9 n + n _i\ 
- sm { J = 2sm [y) COS { 2 J ' ^ 

Writing the right-hand side as a sum of sines, we reproduce equation (125!) with h replaced 

with h = h{l + h 2 /2)- 1 ' 2 . 

The models ( |25l) . (|26l) . (1271) and (1301) constitute our list of new exceptional periodic 

discretisations of the sine-Gordon equation. This list can be generalised and extended in a 

variety of ways. For example, we can replace the sine function in (I2~T1) with tan[(0 n+ i — 9 n )/£] 

or, more generally, with sin[(# n+ i — 6 n )/£] cos p [m(9 n+ i — 6 n )} with arbitrary m and p. Also, 

we can add a sum A n sin 2 [i? n (a; — y)} with arbitrary A n and B n to any of the symmetric 

functions H 2 (x,y). Since we are mainly interested in simple discretisations, we are not 

pursuing these possibilities in our present work. 



V. DISCRETE SINE-GORDON EQUATION WITH EXACT KINK SOLUTIONS 

In this section we consider one more exceptional discretisation of the sine-Gordon equa- 
tion. In addition to admitting an arbitrary centring relative to the lattice, the kinks in this 
model are available in exact explicit form. 

We start with what may seem to be an unrelated map, 

n+ l ~(j) n = h(l- 0„0 n+ l). 

Writing the square of this map as 

(0n+l - 0n) 2 ,2/1 , , \ 

— = h (1 - <p n <pn+l), 

1 - n+ 10„ 

and subtracting its back-iterated copy gives 

4>n+l - 20 n + 0„_i + 4> n {<fil ~ 0n+10n-l) = ~h 2 (f) n (l ~ n+ i0„)(l - 0„0„-l). (32a) 

Equation (I32al) with 

n = COSy, (32b) 

that is, equation 

9n+l 9 n On-l , &n ( 2@n @n+l ^n-l 

cos 2 cos h cos h cos — cos cos COS 



2 2 2 V 2 2 2 

—n cos — 1 — cos — cos 1 — cos — cos , (33) 

2 \ 2 2/V 2 " ' 
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provides an exceptional discretisation of the sine-Gordon equation. Indeed, the continuum 
limit of equation (132aj) is 



) 



(34) 



which is nothing but the stationary sine-Gordon equation 6 XX = sin# written in terms of 
= cos(0/2). 

Equation (132 ap has an exact kink solution 

4> n = tanh(fcn — x^), tanhA; = h, 

where x^ is a translation parameter which can be chosen arbitrarily. Applying the transfor- 
mation (132bj) to <f) n produces an explicit kink solution of the discrete sine-Gordon equation 
(EH}: 

9 n = 4 arctan[exp(fcn — x^)} } tanhA; = h. (35) 

We are not aware of any physical systems represented by Eq. (1331) . This discrete model 
may find its uses, however, in numerical simulations of the continuum sine-Gordon equation. 
Like other exceptional discretisations of the sine-Gordon equation, this model preserves an 
"effective translation invariance" of the continuum equation. The fact that the station- 
ary discrete kinks of the model (133]) are available in exact explicit form is an additional 
computational advantage. 



VI. TRAVELLING KINKS 



In this section we show how the method of one- dimensional maps can be used to construct 
moving kinks. 

The discretisation breaks the Lorentz invariance of the continuum model (T5J) in the same 
way as it breaks its translation symmetry; hence the mobility of the kink becomes a nontrivial 



property in the discrete case. As the kink moves in the 



resonant radiation and decelerates as a result of that 



12 



eier 



14 



s-Nabarro potential, it excites 
20, 



28 



29| | . Surprisingly, some 



discrete models exhibit isolated values of the kink velocity for which the kink can slide, i.e. 



19 



30. 



3l|. 



travel without experiencing radiative friction [18, 

A pertinent question here is whether there are exceptional discretisations where the kink 
can slide with an arbitrary velocity. A discrete nonlinear Schrodinger equation with this 
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property is well known; it is the Ablowitz-Ladik model whose solitons are radiationless 
irrespective of their velocities. On the other hand, no discrete Klein-Gordon equations 
whose kink velocities would all be sliding velocities have been found so far — neither in the 
Frenkel-Kontorova class of models nor among the discrete 4 -theories. 

In this section we construct such a discrete sine-Gordon equation. Its kink solutions are 
given by explicit expressions, and, as will become obvious from these explicit formulas, all 
its kinks travel without emitting radiation. 

We start with a nonstationary equation 

<j> xx - fa - 20^-4 = 0^4, (36) 

which transforms into the sine-Gordon equation 9 XX — 9 tt = sin# by the substitution = 
tan(0/4). Our first observation is that if (j){x,t) is a simultaneous solution of two first-order 
equations 

(37) 



vT 



v 



2 



and 

<Pt = —7^=j4>, (38) 
then it also satisfies Eq. flBU]) . In (J3"7j) and (13"%1) . v is a parameter; — 1 < v < 1. Note that 
Eq. (1361) does not contain v explicitly; hence finding a solution of equations (1371) and ( 1381) for 
all v amounts to finding a one-parameter family of solutions to ( 1361) . 
Next, we discretise Eq. (l3~Tj) according to 



h z 1 — v z 

and divide both sides by the same expression to get 

1 (0n+l - 0n) 2 1 0n0n+l 



(39) 



v (i + C^Xi + K) i - ^ 2 (i + 4>l+i)(i + K) ■ 

We also consider a discrete version of Eq. (158]) : 



(40) 



vT-i; 



2 



(41) 



Subtracting Eq. fjlO]) from its back-iterated copy and replacing <p n with tan(0 n /4), gives 

V Sm { 1 ) = 2 Sm U J C ° S V 4 ) ■ (42) 
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On the other hand, Eq. (T4Tj) yields 



K - 2(fi r , 



hi-Wn+1 



1 + 02 l-v 2 1 + n _i0 n+ i 



where we have used the relation (fi n+1 (fi n _i = <fi n which is straightforward from 
(fi n = tan(6* n /4) in (j4"3l) . we get 



n ( &n+l — On-l 

T cos l i 



1 v 2 



2T^ sm U l<os 



9 n 



(43) 
Letting 

(44) 



Finally, subtracting ff42|) from ff44|) yields a discrete sine-Gordon equation 

^n+l-^n-l\ 1 . f 6 n +l ~ 26 n + 6 n -\\ 1 . / 6 n \ f 9 n +l + ®n-\ 

cos I — = — sin sin — cos 

4 A h 2 \ 4 /2V2/V 4 



(45) 



For any h, this equation has an explicit moving kink solution which is a compatible solution 
of the first-order difference equation ( 1391) and the first-order differential equation (j4"Ii) : 



0„ = 4 arctan 



where is defined by 



2 sinh — 
2 



exp A;n 



VT 



(46) 



(47) 



and v can take any value between —1 and 1. As h —>■ 0, the solution (f46l) -( l47j) tends to the 
travelling kink solution of the continuum sine-Gordon equation, 



9 n — > 4 arctan 



exp 



Using an identity ( 1291) . equation ( 1451) can be cast in the form 



COS I ; Or 



4 . / 9 n +i — 29 n + n -i 

— sm 

h 2 V 4 



6 n+1 + 26 n + n -i i , KM 
sin I : I , (48) 



4 J h? \ 4 J V 4 

where /i = /i(l + h 2 /A)^ 1 ! 2 . Solution to f j48l) is given by the same Eq.( f46l) where k should 
now be defined by 

/ h\ 1 h 

(49) 



sinh I — 



2) VI -v 2 ^A-h 2 
Solution (SSD,(Sn]) exists for any \v\ < 1 and < h < 2. 

We close this section by noting that the stationary limit of Eq.( l48l) coincides with the 
stationary part of the Speight- Ward model [171 ]. 



-i ( @n+i — 9n-i\ -a 4 . / 9 n +i — 26 n + 8 n -i\ . ( n+ i + 26 n + n -x \ 
c,, s ( _ je n = j- 2 sml i 1-sml i ). (50) 
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(Here cos -1 a should be understood as 1/ cos a and not as arccos a.) The simulations of 
Speight and Ward [17J have demonstrated that the motion of the kink in their Eq. (l5UI) is 
accompanied by a much weaker radiation than the kink propagation in a typical nonexcep- 
tional model. Now that we have another time-dependent version of the same stationary 
model, in which the radiation is completely suppressed for all velocities, the low level of 
radiation from the moving Speight- Ward kink can be explained simply by the proximity of 
their equation floUj) to our model fHHj) . 



VII. CONCLUDING REMARKS 



Results of this work can be summarised as follows. 



Using the method of one-dimensional maps, we have derived several exceptional dis- 
cretisations of the sine-Gordon equation involving ratios of trigonometric to linear 
functions: equations ( 1T4"1) . ( floT) . ( fl8i) and ( |20|) . All these exceptional models are new. 
We have also recovered the exceptional system of Kevrekidis, Eq. (112ft . which was orig- 



inally obtained within a different approach 25]. 



We have identified several new purely-trigonometric exceptional discretisations, in 
particular equations ( 1251) . ( |26|) and ( |30i) : 

■a 2 . / 6 n+ i — 29 n + 6 n -\ \ . / 9 n+ i + Q n -\ 
e n = V2 sm ^ j-sm 

I 6n+l — 29 n + 6 n -l \ . I dn+l + dn-l \ ( @n+l ~ @n-l 

6„ = — sin — sin cos 



h 2 V 4 /v 2 yv4 

■a 4 [0 n+ i — 28 n + 9 n -i\ 1 {9 n+ i + 9 n -x\ ( . 9 n+ i . 9 n . 9 n _ 1 
9 n — — sm — - cos sm — — h 2 sin — + sin 



h 2 V 4 / 2 V 4 /V 2 2 2 

We have derived a new discretisation with exact explicit kink solutions, Eq. (l33l) . 

We have constructed a discrete sine-Gordon model which supports kinks travelling 
with arbitrary velocities: 

, ^n+i-^n-i\ 4 f 9 n+ i — 29 n + 9 n _i \ (9 n+ i + 29 n + 9 n _ 1 
"n cos = — sm -sin 



4 / h 2 \ 4 /V 4 



(51) 
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The latter result deserves an additional comment. By analogy with the derivation of the 
model (|5ip . it is not difficult to construct discrete 6 4 theories supporting sliding kinks with 
arbitrary velocities. One such model has the form 

l-0 n+ 10 n _l n+ l - 20 n + n _i 6 n . 
^ n lZT^2 = + y(l-0 n+ i0 n _i). (52) 

This is a time-dependent generalisation of the exceptional stationary 6 4 model derived in 
271].) The moving kink solution to Eq. ( 152ft has the form 

, vt \ 4tanh 2 fc h 2 

tanh kn — 



2y/l - v 2 J 1 + tanh 2 k l-v 2 
Another time-dependent discretisation of the 4 theory with sliding kinks is 



1 - n+ i0 n _i _ (f) n+1 - 2(j) n + n _j (fin 

1-61 h 2 2 



1 _ -Tr(0n+1 + 4>n-l) 



(53) 



(This is a time-dependent generalisation of the exceptional 6 A model identified by Bender 



and Tovbis 32| and Kevrekidis 25].) The sliding kink solution has the form 



i 1 vt \ 2 h 2 

tanh kn — — , , 4 tanh k 



2y/l + tanh 2 A; Vl ~ v 2 J ' ' l-^ 2 ' 

The analogy between equation (15T!) and the Ablowitz-Ladik model is also worth comment- 
ing upon. The Ablowitz-Ladik model is the only discrete nonlinear Schrodinger equation 
whose solitons can slide with any chosen velocity. The absence of the accompanying radia- 
tion is usually explained by the integrability of this equation. A similar behaviour of kinks 
of Eq. (151j) makes one wonder whether the latter equation could also be integrable. We have 
tested the integrability of Eq.( l5T|) numerically, by simulating a collision of a kink and an 
antikink. The scattering was found to be inelastic: the velocities of the kink and antikink 
changed as a result of the collision, and significant amount of radiation was detected. Con- 
sequently, we conclude that equation (jBTj) is not integrable. This example demonstrates 
that, contrary to common belief, the integrability is not a prerequisite for the existence of a 
discrete soliton sliding at an arbitrarily chosen velocity. 

Finally, it is interesting to compare travelling kink solutions of our model ( l5~li) with 
travelling kinks of another modification of the Speight- Ward model proposed by Zakrzewski 



3l|. Zakrzewski's model is different from the Speight- Ward equation (|5T)1) in the presence of 
a factor (1 + aO 2 )^ 1 in front of the left-hand side of (1501). with a = const. The model has an 
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exact solution in the form of a sliding kink; however, similarly to radiationless moving kinks 



in other systems [18|, [19|, |30j, this kink can only slide with one particular velocity which is 
determined by the parameters h and a. Unlike this codimension-1 solution, our sliding kinks 
( |46l) . ( )49i) have codimension in the sense that they can move with an arbitrary velocity 
independent of h. 
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